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Turbulent  Heating  and  Transfer  in  the 
Stratosphere  and  Lower  Mesosphere 


1.  INTRODUCTION 

The  problem  of  determining  turbulence  in  the  upper  atmosphere  up  to  the  turbo¬ 
pause  is  crucial  to  the  understanding  of  this  region.  Not  only  does  turbulence 
transfer  material  and  heat,  but  it  is  also  a  significant  direct  source  of  momentum 
and  heating  to  and  from  the  mean  motions  and  the  thermal  mass  of  the  background 
atmosphere.  These  heating  and  transfer  rates  and  the  mean  height  and  temporal 
behavior  of  turbulence  must  be  ascertained  before  we  will  be  able  to  accurately 
comprehend  the  spatial  and  temporal  distribution  of  species,  temperature,  and 
motions  of  all  forms  (tidal,  planetary,  and  mean)  in  the  mesosphere  and  lower 
thermosphere.  Until  then,  models  are  generated  using  some  known  or  estimated 
forcing  and  damping  functions  to  give  a  good  fit  to  limited  observations. 

In  this  report,  we  attempt  to  define  some  of  the  major  unknowns  of  atmospheric 
turbulence  that  may  be  important  sources  and  sinks  of  local  heating  and  transfer 
of  both  heat  and  momentum.  Preliminary  results  leading  to  this  work  were  reported 
at  the  1979  MAP  Symposium.  There  we  demonstrated  correspondence  between 
analysis  and  some  required  sources  and  sinks  of  heat  and  momentum  in  the  80  to 
90  km  region.  However,  uncertainties  in  the  derivation  of  the  normalized  Reynolds 

(Received  for  publication  22  July  1985) 

1.  Zimmerman,  Z,  P. ,  and  Keneshea,  T.  J.  (1981)  Turbulent  heating  and  transfer 
in  the  stratosphere  and  upper  mesosphere.  Handbook  for  MAP,  2:311-322. 


stress  led  us  to  refine  that  work.  The  principle  question,  we  felt,  was  the  validity 
of  assuming  that  the  ratio  of  the  turbulent  Reynolds  stress  to  the  mean  velocity 
squared,  using  bounding  layer  data,  is  invariant  with  altitude  and  only  a  function  of 
the  local  stability.  We  now  address  the  problem  of  determining  an  altitude-invariant, 
normalized  and  non-dimensional  turbulent  parameter. 


Z/L 


Figure  1.  The  Schematic  Representation  of  the  Dimensionless 
Energy  Budget  Under  Unstable  Conditions  (after  Figure  7, 
Wyngard  and  Cote,  197 1^) 


We  approach  the  problem  of  describing  atmospheric  turbulence  by  considering 
the  turbulent  energy  budget.  We  shall  show  that  in  the  boundary  layer,  for  both 
statically  stable  and  unstable  atmospheric  regions,  the  concept  of  volumetric  rate 
of  energy  balance  of  Reynolds  stress  (<u'  heat  flow  <w'  d'>),  and 

dissipation  (e)  provides  a  fairly  good  description  of  the  local  turbulent  region,  where 
the  sum  of  Reynolds'  stress,  heat  flux,  and  dissipation  is  zero. 
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is  the  turbulent  Reynolds'  stress, 
is  the  turbulent  heat  flux, 

is  the  rate  of  dissipation  of  turbulent  kinetic  energy  by 
kinematic  viscosity, 

is  the  local  vertical  gradient  of  the  horizontal  wind  (U), 


is  the  mean  potential  temperature  =  T  (-r^) 


denote  ensemble  averages. 


For  this  study,  we  have  selected  what  we  consider  a  more  logical  function  for 
extrapolation  from  the  boundary  layer  to  the  mesosphere  than  those  previously  used.  ^ 
This  function  is  the  non-dimensional  Reynolds' stres  spectral  length  scale  normalized 
to  the  wind  length  scale.  By  demonstrating  the  invariance  with  altitude  of  this  non- 
dimensional  ratio,  we  can  apply  it  to  other  bodies  of  wind  and  temperature  data  to 
determine  the  turbulent  parameters.  Using  this  analysis  we  present  probable  values 
of  the  turbulent  parameters  and  their  altitude  variations  from  the  stratosphere  and 
into  the  mesosphere. 


2.  THEORY 

First  let  us  examine  the  point  we  raised  in  the  introduction  about  the  assumption 

of  rate  of  energy  balance  between  turbulent  source  and  heat  loss,  or  production  and 

2 

dissipation.  Wyngaard  and  Cote  demonstrate  that,  for  boundary  layer  data  and  a 

2.  Wyngaard,  J.  C.  ,  and  Cote,  O.  R.  (1971)  The  budget  of  turbulent  kinetic  energy 
and  temperature  variance  in  the  atmospheric  surface  layer,  J.  Atmos.  Sci. , 
28:190-201. 


positive  Richardson  number  this  rationale  is  well  established.  However,  for 

negative  Rj,,  their  values  show  that  generally  the  balance  is  nearly  maintained  with 
the  energy  flux  terms  (Figure  2).  This  is  also  shown  by  Caughey  and  Wyngaard 
in  a  separate  work  (Figure  3).  Thus,  the  general  assumption  of  energy  rate  balance 
that  we  have  made,  is  fairly  well  substantiated  in  lower  atmospheric  studies. 
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Figure  2.  The  Same  as  Figure  1  Except  that  the  Imbalance 
Has  the  Turbulence  Transport  Term  Subtracted  From  it. 
This  suggests  that  even  for  negative  stability,  dissipation 
is  approximately  balanced  by  buoyant  production  plus  shear 
production 


We  now  consider  the  spectral  relations  for  the  Reynolds' stress  and  heat  flux. 

4  5 

Following  Wyngaard  and  Cote,  and  Kaimal  et  al  the  Reynolds'  stress  is 
expressed  as 


3.  Caughey,  S.  J, ,  and  Wyngaard,  J.  C.  (1979)  Some  aspects  of  turbulence  structure 

through  the  depth  of  the  boundary  layer.  Quart.  J.  R.  Met.  Soc. ,  1  (35:811-827. 

4.  Wyngaard,  J.  C. ,  and  Cote,  O,  R.  (1972)  Cospectral  similarity  in  the  atmos¬ 

pheric  surface  layer.  Quart.  J.  R.  Met.  Soc. ,  pp.  590-603. 

5.  Kaimal,  J.  C. ,  Wyngaard,  J.  C. ,  Izumi,  Y. ,  and  Cote,  O.  R.  (1972)  Spectral 

characteristics  of  surface  layer  turbulence.  Quart.  J.  R.  Met.  Soc.  , 
pp.  563-589. 
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Figure  3.  The  Same  as  Figure  2  but  for 
the  Data  of  Caughey  and  Wyngaard.  ^ 
Components  in  the  energy  equation  non- 
dimensionalized  with  the  mixed  layer 
similarity  scale  vr*  and  z^,  plotted  against 
z/L.  The  full  lines  are  simply  drawn 
through  the  data  points. 
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and  the  heat  flux  as 


/  ,0,-^  ,1/3  00  r, -7/3 

<w'0'>=  -  Fi  £  '  57  7  k 


and  where  the  integration  is  over  the  inertial  range  between  wavenumber  k^,  the 
energy  containing  eddy,  and  wavenumber  k^,  the  dissipative  scale.  Cj  and  are 


I 


spectral  coefficients  that  are  functions  not  only  of  Z/L  (the  ratio  of  altitude  to  the 
Monin -Obunkhoff  length  scale),  but  also  of  the  local  volumetric  Richardson  numbers 
(Kaimal,  1973).  The  integral  is  taken  only  over  the  inertial  spectral  function,  since 
shear,  dissipation,  and  temperature  gradient  are  assumed  independent  of  wavenumber. 

Now  let  us  consider  the  simple  solution  for  the  heat  flux  by  integrating  over  the 
spectrum  from  k^  (the  outer  scale  of  the  turbulent  spectrum)  to  <».  Then  Eq.  (2) 
becomes 

^  3  1/385,  -4/3 

<w'5->  =  -  ^  k^  '  .  (3) 


Solving  for  l/k  )  with  the  assumption  of  rate  of  energy  balance  and  the  norma¬ 
lization  to  the  vertical  length  scale  of  the  mean  horizontal  wind  {I  =  U/|^)  we 
arrive  at 


<w'  d'> 
~l/2  , 


„  1  ,e  85/'^  -2 


^  1 


Where  is  the  gradient  Richardson  number  ^ 
Similarly  for  the  Reynolds'  stress,  we  may  show 


■fuA 


87 \2' 


(1  -  Rj)^'^  (7)' 


4  5 

Wyngaard  and  Cote  and  Kaimal  et  al  have  shown  that  the  heat  flux  spectral 
coefficient  Vj  and  the  stress  spectral  coefficient  Cj  are  functions  of  Z/L  from  4.  to 
26  m  altitude.  However,  in  the  upper  atmosphere  the  Monin-Obukhoff  lengths  or 
the  boundary  layer  height  cannot  be  used  as  scaling  factors.  Thus  we  must  examine 
other  concepts.  Our  choice  was  to  find  a  volumetric  scaling  factor  and  then  to 
derive  the  volumetric  relations  of  y  and  Cj,  as  functions  of  Rj.  Then  by  com¬ 
paring  these  determined  parameters  to  the  above  authors  measurements,  we  can 
determine  the  unknown  constants. 

To  solve  for  these  coefficients,  we  use  Eq.  (3)  and  its  equivalent  for  the 

Reynolds'  stress,  and  substitute  the  fluctuating  buoyancy  conditions  for  the  outer 

scale  I  ,  where 
o 

-1  r 

L'L  -=1  ''4  "" 
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where  u)g  is  the  Brunt-Vaissala  frequency.  ^  Wg  =  ^  ^  J  , 


Again  using  rate  of  energy  balance  to  solve  for  6,  Vj  and  it  can  be  shown 
(Appendix  A)  that 


Rr 


y  1  =  - 


3  C, 


^  j  .  -2.  5  <  Z/L  <  -0.  3 


(7) 


and 


^1  = 


3  C. 


1  -  R. 


-2.  5  <  Z/L  <  -0.  3 


(8) 


for  statically  unstable  conditions.  The  constant  Cj  is  determined  by  normalization 
to  Figure  2  of  VVyngaard  and  Cote^  (our  Figure  4,  open  circles).  The  results 


(F  igure  4)  show  that  for  more  negative  R^,  =  1.  66  using  Eq.  (7)  and  C2  -  !•  2 

using  Kq.  (8). 
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Figure  4,  The  Vertical  Heat  Flux  Spectral  Coefficient 
a  j  as  a  Function  of  Z/L.  The  open  circles  are  the 
data  of  Wyngaard  and  Cote,  ^  and  the  solid  circles  from 
the  relations  Eqs.  (A5)  and  (A7> 


For  the  near  neutral  and  stable  region,  say  -0.  3  <  <  0.  3,  rate  of  energy 

balance  has  the  approximate  relation 
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<u'  w'> 


au 

5F 


(9) 


and  using  Eq.  (A3)  (Appendix  A,  Eq.  (A7)),  we  observe  that 


4 


(10) 


By  normalization  to  the  data  Cg  =  0.  65,  and  using  a  similar  approach  for 
<u'  w'>,  we  see  that  [Appendix  A,  Eq.  (A  10)) 


C4  =  0.  35  .  (11) 

is  displayed  in  Figure  5  for  -2.  5<Z/L<  3  which  is  approximately  equivalent 
to  -2.  5  <  Rj  <  0,  3. 
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Figure  5.  The  Same  as  Figure  4  but  for  the  Reynolds' 
Stress  Spectral  Coefficient  C, 


As  observed  there  is  quite  good  agreement  between  y  ,  measured  by  Wyngaard 

4  * 

and  Cote  and  the  derived  relation  Eq.  (7)  for  Rj  between  these  limits.  For  values 

of  Rj  between  -0.  3  and  +0.  2  (-0.  3  <  Z/L,  +2,  5)  the  above  expression  does  not  fit 

the  observations.  However,  the  empirical  relation  Eq.  (11)  describes  these  results 

fairly  accurately.  The  comparison  between  Eq.  (8)  and  the  data  points  determined 

0 

from  the  Kansas  experiment  are  rather  scattered,  with  both  analyses  showing  sig¬ 
nificant  fluctuations  around  the  Kaimal  et  al  empirical  curve.  Given  the  relation 
for  y  .  and  ?  ,  we  now  examine  the  non-dimensional  length  scales  (4)  and  (5) 

11  .j 

using  both  the  Kansas  and  the  Minnesota  data  bases.  Figures  6a  and  6b  are  the 
representation  of  Eq.  (5)  for  the  Kansas  data  and  the  Minnesota  data  respectively. 

We  see  that  the  large  scatter  exhibited  in  the  Cj  coefficient  is  not  so  apparent  in 
this  non-dimensional  form  of  the  Reynolds'  stress  for  both  the  unstable,  neutral 
and  stable  regions.  The  values  displayed  in  Figure  6b  show  much  greater  scatter 
but  they  follow  quite  closely  the  Kansas  results.  Figures  7a  and  7b  show  the  values 
from  Eq.  (4)  for  the  same  conditions  as  described  for  Figures  6a  and  6b.  Again 
we  observe  quite  reasonable  results  utilizing  both  bodies  of  data.  Thus,  we  have 
demonstrated  the  existence  of  a  non-dimensional  description  of  the  major  compo¬ 
nents  of  the  turbulence  rate  of  energy  balance  equation  from  4  meters  to  the  top 
of  the  boundary  layer.  The  observation  that  the  Minnesota  results  follow  quite  well 
and  bracket  the  Kansas  results  suggests  that  these  non-dimensional  scales  are  a 
viable  way  of  describing  turbulence  through  the  boundary  layer.  It  is  definitely 
not  a  complete  description,  since  we  assume  the  pressure  term  and  the  energy  flux 
term  completely  balance.  But  for  our  purpose,  the  description  of  turbulence  through 
the  troposphere  to  the  mesosphere,  this  assumption  may  be  fairly  valid.  To  further 
ascertain  the  validity  of  our  results,  we  examine  other  data  samples  obtained  by 
aircraft  up  to  20  km  (Figure  8).  As  observed,  the  non-dimensional  heat  flux 

data  points  fall  quite  nicely  around  the  calibration  curve  (the  Kansas  boundary  layer 
date)  to  approximately  20  km.  Given  this  normalization,  we  feel  that  the  extrapola¬ 
tion  of  these  relations  to  the  approximate  turbopause  (say  100  km)  is  reasonable. 

(Due  to  the  large  number  references  cited  above,  they  will  not  be  listed  here. 

See  References,  page  25.  ) 


Figure  6.  The  Representation  of  the  Non-dimensional  Reynolds’  Stress  [Eq,  (5)j 
Using  the  Boundary  Layer  Data  From:  (a)  the  Kansas  Experiment  (Izumi®)  and 
(b)  the  Minnesota  Experiment  (Izumi  and  Caughey^)  Respectively.  As  observed, 
there  is  no  systematic  altitude  dependence  observed  in  the  analyzed  data,  and 
both  sets  follow  quite  similarly  to  each  other,  for  both  positive  and  negative 
stability.  The  Kansas  results  show  less  variability  than  that  of  the  Minnesota 
analysis.  We  think  that  this  is  due  to  the  fine  measuring  grid  used  in  the 
Kansas  experiment.  This  would  result  in  better  determination  of  temperature 
and  wind  gradients  and  thus  more  accurate  values  of  those  quantities  used  in 
these  analyses 


Kansas 

■  5.  66  m. 

+  1 1.  3  1  m, 
O  22.  63  m. 


Minnesota 

□  1219  m. 
A  9  14  m. 
O  610  m. 
H  457  m. 
^  305  m. 

X  150  m. 


3.  PREVIOUS  ANALYSES  AND  CONCEPTS 


We  now  analyze  for  the  turbulent  transfer  coefficients  and  address  here  two 

problem  area.  First  is  the  question  of  the  critical  flux  Richardson  number  (R^  ) 

and  what  is  its  value  when  the  ratio  of  heat  energy  transfer  by  turbulence  to  direct 

turbulent  heating  is  unity.  This  question  has  been  discussed  in  the  literature  by 
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many,  including  work  by  Hunten,  and  Izakov.  We  shall  demonstrate  here  that, 
contrary  to  previous  works  that  infer  a  constant  number  near  unity,  this  ratio  has, 
within  limits,  an  infinity  of  numbers. 

The  second  problem  area  concerns  the  effects  that  turbulence  has  on  the  mean 
thermodynamic  systems.  These  are  the  rate  at  which  turbulence  degrades  or 
enhances  the  mean  thermal  field,  and  the  production  and  loss  of  momentum  by  tur¬ 
bulence  to  and  from  mass  motion.  In  the  mesosphere  these  problems  have  been 
addressed  by  Ebel^^'  and  by  Chandra.  Ebel  analyzed  Groves' model 

atmospheres  for  the  heating  and  momentum  supply  and/or  loss  that  must  come  from 
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external  sources.  Chandra  following  Izakov's  analysis,  theoretically  examined 
the  conditions  in  the  mesosphere  and  upper  stratosphere  for  local  turbulent  heating 
by  dissipation  and  cooling  by  turbulent  transfer,  in  a  steady  state  energy  balance 
model.  He  suggests  that,  for  a  given  flux  Richardson  number  (R^  =  0.  23  =  R^  ), 
the  ratio  of  cooling  by  turbulent  transfer  to  heating  by  dissipation  equals  unity? 

For  Rj  <  0.  23  the  ratio  is  less  than  unity  and  for  R^  >  0.  23  the  ratio  is  greater  than 
unity.  Chandra  then  utilizes  a  value  of  =  0,  375  to  calculate  the  temperature 
structure,  and  heat  production  and  loss,  of  the  upper  mesosphere  and  lower 
thermosphere.  We  shall  examine,  in  the  light  of  the  work  presented  here,  Ebel's 
and  Chandra's  results  as  well  as  Chandra's  assumption.  We  shall  also  present 
determined  values  of  turbulent  heat  and  momentum  production,  loss  and  diffusivities. 

4.  THEORY  OF  THE  MEAN  AND  TURBULENT  HEAT  TRANSFER 

We  will  now  show  how  these  processes  are  determined  in’ the  atmosphere.  Let 
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us  start  with  the  heat  energy  equation,  Eq.  (12)  Chapman  and  Cowling,  Banks  and 
Kockarts,  and  Izakov. 


(Due  to  the  number  of  references  cited  above,  the  will  not  be  listed  here. 
See  References,  page  25.  ) 
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is  the  total  heat  flux  vector, 

is  the  mass  density, 

is  the  specific  heat  at  constant  volume, 

is  the  specific  heat  at  constant  pressure, 

is  temperature  (K) 
is  pressure 

are  the  chemical,  turbulent,  and  radiative  heat  source  and  loss, 

is  the  adiabatic  lapse  rate, 
is  the  acceleration  of  gravity, 
is  the  mean  mass  velocity 

is  the  mass  of  ith  species, 

is  its  average  velocity,  and 

is  the  turbulent  momentum  diffusivity  tensor. 


The  turbulent  component  can  be  shown  (Appendix  B)  to  be  equal  to  the  divergence  of 
the  vertical  heat  flux  term  <w'  T'>  which  is  the  correlation  of  the  vertical  velocity 
and  temperature  fluctuations. 
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where 


<w'  T>  is  the  turbulent  kinetic  temperature  heat  flux, 
w'  is  the  vertical  component  of  turbulent  velocity,  and 

T'  is  the  turbulent  temperature  fluctuation. 
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Similarly,  from  the  equation  of  motion 
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Where  we  ignore  the  nonlinear  and  Coriolis  terms,  we  see  that  the  turbulent 
vorticity  expression  reduces  to  the  divergence  of  the  turbulent  Reynolds  flux. 
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Here  the  positive  and  negative  values  of  the  turbulent  momentum  transfer  are 
directed  along  the  wind  vector,  and 

C  is  the  mean  wind  vector, 
o 

p/p  is  the  kinematic  viscosity,  and 
p  is  the  molecular  viscosity. 

Finally,  by  going  to  the  steady  state  turbulence  energy  rate^q.  (16)  and  using 
the  expression  for  the  flux  Richardson  number  (R^) 


-  <u'  w'>  ^  <w'  T'>  =  € 


we  can  express  the  ratio  of  the  heat  transfer  divergence  to  turbulent  dissipation 
as  (Appendix  C) 


Ratio  = 


f  j  ^  “  <w'T'  > 


where  and  are  the  density  and  heat  flux  scale  heights.  Results  from 

this  analysis  are  shown  in  Figure  9. 

We  now  have  relations  for  determining  the  local  production  and  loss  of  heat, 
momentum  and  the  rate  of  turbulent  dissipation  from  rate  of  energy  balance.  We 
may  determine  the  diffusivities  also,  since  we  have  knowledge  of  the  mean  tempera¬ 
ture  and  wind  fields.  These  last  parameters  are  simply  related  by 


<w'T->=  + 


<u'w'>  =  -K^ 
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(19) 


where  K_  and  K  are  the  turbulent  diffusion  coefficients  of  heat  and  momentum 
T  m 

respectively. 


ROTIO  vs  AF  MALLOPS  ISLRNO 


Figure  9.  The  Ratio  of  the  Turbulent  Heat  Flux  Divergence 
to  the  Rate  of  Turbulent  Dissipation  (Ratio)  as  a  Function  of 
Rf.  Also  shown  are  theoretical  curves  for  different  values  of 
Ho/H<w'  0<  >  (the  ratio  of  the  density  scale  height  to  the 
turbulent  heat  flux  scale  heights).  The  data  to  the  left  of  the 
positive  vertical  line,  at  the  ratio  equal  to  unity,  marks  the 
transition  area  where  cooling  by  turbulent  transfer  dominates 
over  the  heating  by  dissipation.  The  remainder  of  the  area 
denotes  effective  heating.  The  right  hand  side  denotes 
local  heating  both  turbulent  transfer  and  dissipation 


5.  DATA 


The  data  base  used  to  determine  the  turbulence  coefficients  and  their  winter- 
summer  altitude  and  latitude  variability  are  the  measurements  from  meteorological 
rocketnetwork  (MRN)  covering  the  period  1969  to  1978.  The  data  consists  of  wind, 
temperature  and  some  density  measurements,  generally  at  height  intervals  of  1  km. 
The  data  were  splined,  and  the  data  points  were  interpolated  at  500-m  intervals, 
at  which  point  the  derivatives  of  wind  and  temperature  were  determined.  Figure  10 
is  a  sample  temperature  profile  showing  the  measurements,  the  500-m  spline 
values  and  mid-point  average  values. 
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Figure  10.  An  Example  of  the  Temperature  Measurements  from  18  to  70  km. 
The  ordinate  and  abclssa  of  the  measurement  are  divided  into  10  divisions. 
Superposed  on  the  graphs  are  the  mid-point  average  (+)  and  the  mid-point 
spline  (x) 


These  data  are  from  a  group  of  observational  sites  covering  the  latitude  band 
76.  6'  N  to  9.  3“  N  and  the  narrow  longitude  band  from  75.  5°  W  to  80.  5®  W.  In  this 
longitudinal  band,  extending  from  the  northern  most  point  to  as  far  south  as  a  sub 
stantial  data  base  allows,  the  variability  of  the  turbulence  dynamics  as  a  function 
of  northern  latitude  winter  and  northern  latitude  summer  is  determined. 


m 


6.  RESULTS 


Based  on  the  rocket  data  at  Fort  Sherman  (9.  3°  N  Lat  and  80°  W  Long), 

Figure  11  shows  the  statistics  of  the  Richardson  number  for  bin  widths  of  0.  0625. 
The  important  feature  here  is  that  the  most  probable  R/s  are  in  the  region  0  to 
0.  187  as  contrasted  to  that  of  the  stratosphere  wherein  the  most  probable  value  is 
considered  to  be  0.  25.  This  significant  difference  reflects  strongly  the  difference 
between  stratospheric  and  these  lower  mesospheric  microdynamics.  This  would 
specify  that  for  a  stably  stratified  turbulent  region,  the  deduced  turbulence  param¬ 
eters  (K  and/or  £  )  determined  from  the  rate  of  energy  balance,  would  be  greatly 

changed  if  the  value  of  R.  used  was  0.  05  rather  than  0.  25.  This  has  further  rele- 
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vance  with  regard  to  the  energy  balance  calculations.  Chandra  assumed,  as  a 

first  guess,  that  the  effective  value  of  R^  was  ~  0.  375  for  the  entire  stratosphere 

and  mesosphere.  This  would  be  a  highly  stable  atmosphere,  and,  he  points  out, 

would  insure  that  in  his  calculations  the  divergence  of  the  turbulent  heat  flux  term 

would  always  be  a  positive  cooling  mechanism  over  the  entire  mesosphere.  A  shift 

to  a  lower  R^,  say  0.  05  in  the  7  0  to  90  km  region  of  the  mesosphere  would  introduce 

significantly  more  heating  into  this  region  by  markedly  increasing  the  turbulent 

dissipation  rate  for  an  adopted  diffusivity.  However  in  these  data,  the  other  sites 

show  the  most  probable  R^^  to  be  larger  than  0.  25  and,  under  these  conditions 

(R.  s  0,  25),  Chandra's  analysis  is  more  appropriate. 

The  values  of  the  turbulent  parameters  determined  from  these  MRS  measure¬ 
ments,  and  within  the  longitudinal  restraints  we  have  imposed,  are  presented  in 
Figures  12,  13  and  14.  The  specific  locations  of  these  measurements  are  Thule, 
Greenland  (70.  6°  N,  68.  8°  VV),  Wallops  Island,  Virginia  (37.  8°  N,  75.  5°  W), 

Cape  Kennedy,  Florida  (28.  5°  N,  80.  5°  W)  and  Fort  Sherman,  Canal  Zone  (9.  3°  N 
and  80°  W). 

The  criteria  we  used  that  limited  the  results  to  the  above  sites  was  that  the 
number  of  measured  observations  must  be  sufficient  enough  to  supply  roughly 
contiguous  results  within  the  measured  altitude  range.  This  enabled  us  to  compile 
a  meaningful  set  of  measurements  to  specify  the  effective  values  of  the  turbulent 
parameters.  The  parameters  analyzed  for  are  the  heat  diffusivity  (K„),  the 

Q  rp  rl 

momentum  transfer  (9u/9t)  and  total  heating  cooling  (^  +  e). 
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Figure  11.  Frequency  of  Occurrences  of  Richardson  Number  for 
Fort  Sherman 


The  values  of  each  of  these  turbulent  parameters,  for  all  sites,  are  shown  in 
a  sequential  arrangement  along  a  latitudinal  path  (Figures  12,  13  and  14)  from 
Thule,  Greenland  to  Fort  Sherman,  Canal  Zone.  The  winter  periods  cover  the 
six-month  period  around  the  winter  solstice.  Some  of  these  parameters  follow  a 
remarkably  sequential  profile  as  a  function  of  altitude  considering  the  nine-year 
time  period  of  measurement  and  the  logic  that  we  are  examining  the  fluctuating 
components  of  wind  and  temperature.  A  case  in  point  is  the  turbulent  diffusivity 
(Kjj)  measured  at  the  Canal  Zone  9.  3°N  Latitude.  As  observed  (Figure  12)  both 
the  winter  and  summer  periods  show  small  variation  from  a  profile  that  is  a 
function  of  altitude  only.  As  observed  (particularly  in  Du/Dt)  there  does  not 
appear  to  be  a  balance  of  energies.  This  is  particularly  so  because  each  parameter 
displayed  is  over  a  band  of  amplitude  that  would  have  some  effect  upon  the  total 
dynamics  of  the  atmospheric  region.  If  the  amplitude  of  any  parameter  fell  below 
the  lowest  limit  of  the  band,  then  it  would  not  be  displayed  in  the  plot.  This  is  the 
reason  for  the  extreme  paucity  of  data  in  the  plots  of  9u/9t  in  the  northern  latitude 
(Thule)  as  compared  to  those  of  the  "Canal  Zone". 
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Zimmerman  and  Keneshea,  using  the  same  analytic  procedure  in  their 
examination  of  the  rocket  grenade  data,  noted  that  the  northern  latitude  (Poker  Flat, 
Latitude  7  0°  N)  winters  are  much  more  energetic  than  the  summers,  and  show  a 
higher  occurrence  rate  of  turbulence.  Generally,  these  findings  are  substantiated 
in  these  analyses,  but  as  observed,  this  trend  is  not  completely  substantiated  when 
we  examine  other  latitudes.  The  Wallops  Island  results  of  do  indeed  show  this 
winter-summer  trend  as  does  Cape  Kennedy,  but  the  Fort  Sherman  (Canal  Zone) 
results  indicate  a  reversal  of  these  dynamics,  where  the  summer  turbulence  is 
more  intense  than  that  of  the  winter. 

When  considering  the  diffusivities  (Figure  12)  we  observe  the  intense 

variability  of  the  northern  latitude  (Thule)  winter,  and  the  more  uniform  behavior 

of  the  summer  season.  This  is  similar  to  the  results  of  the  rocket  grenade  data 
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(Zimmerman  and  Keneshea  ’  )  but,  because  of  the  finer  scale  of  the  measurement 

grid,  we  have  turbulence  measurements  to  lower  altitudes  for  both  winter  and 
summer.  This  is  particularly  so  for  the  lower  latitudes  (Cape  Kennedy  and  the 
Canal  Zone)  where  the  data  shows  an  almost  contiguous  profile  dov  n  to  approximately 

~35  km.  The  Wallops  data  near  30  km  shows  a  ledge  of  amplitude  diffusivity.  This 
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phenomenon  was  also  noted  by  Zimmerman  and  Keneshea.  These  large  amplitude 
energetics  are  also  noted  in  the  turbulent  wind  acceleration  and  total  heat  produc¬ 
tion  and  loss  mechanisms.  Moreover,  at  Cape  Kennedy,  we  also  observe  sig¬ 
nificant  winter  heating  at  the  lower  altitudes,  and  again,  in  the  summer.  There  is  a 
suggestion  in  these  energies  that  they  must  be  related  to  a  mid-latitude  source 
constrained  in  the  upper  stratosphere,  perhaps  the  "jet  stream"  that  migrates 
around  the  mid  latitudes.  Further  analysis  of  more  extensive  sites  should  deter¬ 
mine  the  source  of  these  energies. 

The  turbulence  wind  acceleration  is  generally  balanced  positive  and  negative 

except  for  the  winter  at  Cape  Kennedy,  where  turulence  provides  an  effective  drag 

on  the  wind  system,  and  at  Fort  Sherman  in  the  summer,  where  there  is  a  large 
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net  acceleration  to  the  mean  winds  to  the  value  of  10  m/sec  .  The  total  turbulence 
heating  (Figure  14)  shows  a  general  balance  of  heat  production  and  loss  in  the  mid 
and  northern  latitudes,  but  at  Cape  Kennedy  and  Fort  Sherman  for  the  winter,  and 
around  60  to  70  km,  there  appears  to  be  a  significant  loss  of  heat.  In  the  summer 
however,  at  Fort  Sherman,  turbulence  is  depositing  significant  heat  into  the 
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mesosphere  to  the  value  ~  3  X  10  K/sec  or  ~  300  K/day. 

19.  Zimmerman,  S.  P. ,  and  Keneshea,  T.  J.  ( 1985)  Turbulent  heating  and  transfer 
in  the  stratosphere  and  upper  mesosphere,  J.  Atmos.  Terr.  Phys. , 
to  be  published. 


This  initial  analysis  of  the  MRN  data  for  turbulence,  its  amplitudes  and  sign, 
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substantiates  the  results  of  Zimmerman  and  Keneshea  for  the  Wallops  Island 
data.  The  results  from  the  other  regions  along  this  longitude  band  show  that 
mesospheric  and  upper  stratospheric  turbulence  has  sufficient  energies  to  significantly 
modulate  mean  motions  and  species  distributions. 
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Appendix  A 


The  calculations  in  the  body  of  this  article  are  based  strongly  upon  the  values 

of  the  spectral  coefficients  y  ^  and  as  functions  of  the  Richardson  numbers 

and  R..  The  experimental  variation  of  y  and  C,  are  demonstrated  by  Wyngaard 
i  ^  1  1 

and  Cote  as  functions  of  the  ratio  of  the  height  of  measurement  (Z)  to  the  Monin- 

Obukhoff  scale  (L).  However,  for  the  atmospheric  region  above  the  boundary  layer, 
this  scaling  factor  no  longer  holds  and  we  must  utilize  the  volumetric  Richardson 
numbers  in  the  expressions  of  these  spectral  coefficients.  Since,  as  will  be  shown, 
no  single  relation  yet  explains  the  measurements  for  all  R^,  we  must  utilize  differ¬ 
ent  derived  expressions  for  these  separate,  but  contiguous,  regions  of  Rj. 

To  derive  these  analytic  forms  of  y  ^  and  we  start  with  Eq.  (3), 


<w'e'>  =  -  i  y  ^ 


1/3  9  0  ,  -4/3 
‘  TZ  o 


Now  using  the  outer  scale  relation  Eq.  (6), 


k  C 

o  ~  ‘-I 


This  is  a  fluctuating  condition,  where  is  the  square  of  the  Brunt  Vaissalla 

0  d  ^  ^ 

frequency  (■  ^  then  becomes 


<w'0'>  =  -  C  y  ^ 


Utilizing  rate  of  energy  balance  (<-  u'  w'>  ^  <w'  9  ’>  =  e)  to  solve  for  e,  and 

the  relation, 

■^  <  w'  0  '> 


we  solve  for  y  j,  and  arrive  at 


^  1”  3i 


which  is  Eq.  (7)  in  the  main  text. 


As  we  observe  (Figure  4)  this  relation  normalized  to  the  peak  amplitude  is 
only  descriptive  of  the  region  -7  <  Rj.<  -0.4,  or  -2.  5  <  ^  <  --0.3.  For 
Rj.  >  -0.4,  there  are  large  deviations  from  the  measured  values. 


To  derive  an  analytic  expression  for  this  region  (-0,4  <  R- <  1),  we  again  start 

^  5 

with  Eq.  (Al),  utilize  the  outer  length  scale  condition  (Kaimel  et  al  ),  and  arrive 


and  observe,  using  rate  of  energy  balance,  that 


^  irj  • 


This  relation  Eq.  (A7)does  indeed  represent  y  j,  quite  nicely  for  the  region 


-0.  4  <  Rj,  <  0.  2  |Eq.  (A2)]. 


In  a  manner  similar  to  the  heat  flux,  we  utilize  the  equation  of  the  spectral 
Reynolds'  Stress, 


/  ,  3  ,  1/3  9  U  ,  -4/3 

<u'w'>=-4?i€  g-^k^ 


Substituting  the  fluctuating  condition  [Eq.  (A2)]  for  k^,  and  using  rate  of  energy 


balance,  we  arrive  at 


■^1  T(^  (  1  -  Rj 


which  is  Eq.  (A8)  in  the  text. 


I 


S' 
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Figure  5  shows  this  [Eq.  (A4)]  is  a  fairly  good  description  of  ?  j  for  large 
negative  values  of  Rj,  with  large  deviations  as  R^  -*  0,  Using  Eq.  (A8)  and  condition 
Eq.  (A6),  we  observe  that 


for  the  condition  -0.  3  <  R^  <  1.  The  results  (Figure  5)  show  this  relation  to  be  a 
fair  representation  of  the  measurements,  which  themselves  are  not  in  as  concise 
a  grouping  as  those  of  y  j^. 

These  analytic  functions  are  used  to  supply  the  values  of  y  ^  and  ?  ^  after  the 
local  Richardson  number  has  been  determined  from  the  wind  and  temperature  data. 


a 


Appendix  B 


From  Eq.  (12)  we  consider  only  the  turbulent  terms,  and  ignore  dissipation 
(t ).  Thus,  in  the  vertical  dimension  we  have. 


a  T  _  1 _  a 

TF  "  pO  ST: 

V 


Setting  the  ratio  of  the  specific  heats  equal  to  y  and  using  the  identity 


<w'  T'>  =  -  I  FZ  ^ 


we  have 


a  T 


(p  M 


p  <  w'  T'> 


Then  simply  differentiating,  we  arrive  at 


^  -  y  <  w'  T'>  - — <  w'  T’>  +  -  ^ 
a  t  '  <  w'  T'>  a  Z  0  aZ 


which  is  Eq.  (13)  of  the  text. 


Appendix  C 


The  ratio  of  the  divergence  of  the  turbulent  heat  flux  (<  w'  T'>)  to  the  rate  of 
turbulent  dissipation  has  been  expressed  in  the  text,  (Eq.  (17)].  We  now  consider 
the  derivation  of  this  relation. 

The  ratio  (K)  is  given  by. 


and  by  assuming  rate  of  energy  balance,  pe  may  be  expressed  as 


<  w'  T'>  . 


Expanding  the  numerator  of  Eq.  (Cll  and  rearranging,  we  have  for  the  ratio 


Ij) /  Hj.  <  w'  T'>  5^ 


<  w'  T'>  + 


1  8p 
I)  9Z 


.Now  Cp  the  specific  heat  at  constant  pressure,  for  a  diatomic  molecule  is  equal 

to  Y  Where  k  is  Boltzmann's  constant  and  m  is  the  mass  of  the  molecule. 

Thus  lEq.  (3']  becomes 


y  i.j  ■=>  *<w'  T'>  “p  j 


<w'  T'> 


where  is  the  heat  flux  scale  length,  defined  as  '  T>  “  8 - 

and  Hp  is  the  density  scale  height,  defined  as  •  Since  is  the 

atmospheric  pressure  scale  length  Hp,  we  then  have 


R-1  f_X 
^  2  I  1  -  R, 


<w'  T'> 


We  now  assume  that  Hp  =  H^,  which  specifies  that  Ht  (the  temperature  scale 
length)  »  Hp.  According  to  the  U.  S.  Standard  Atmosphere  (1976)  this  is  a  fairly 
valid  assumption,  throughout  the  atmosphere.  The  ratio  then  becomes 


Vw'  T'> 


This  relation  [Eq.  (C6)],  for  values  of  (Hp/H<w'T'>)  equal  to  -100,  -10,  -2,  0,  2 
10,  and  100,  are  included  in  Figure  9  as  the  solid  curved  lines. 
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